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SUMMARY 


Studies of various arrangements of vings and bodies designed to 
provide favorable wave interference at supersonic speeds lead to the 
problem of determining the minimum possible value of the wave resist- 
ance obtainable by any disposition of the elements of an aircraft within 
a definitely prescribed region. Under the assumptions that the total 
lift and the total volume of the aircraft are given^ conditions that 
must be satisfied if the drag is to be a minimum are found. For arbi- 
trary regions the minimum value can be estimated by a simple formula 
giving a lower bound. 


DISCUSSION 


In 1935 Busemann (ref. l) showed that the wave drag of two airfoils 
could be canceled by reflection. Later Ferrari (ref, 2) showed that the 
drag of a body of revolution could be canceled by the addition of a ring 
airfoil to catch the wave from the nose and reflect it back to the tail. 
Even if the investigation is limited to such completely self-contained 
wave systems^ these examples are only two of an infinite number of 
possibilities . 

The examples in which the wave cancellation is complete are, however, 
limited to systems in which the net lift and lateral force are zero. 
Nevertheless, examples cited by Ferri (ref. 3 ) Graham (ref. 4) show 

that the wave drag associated with the lift can be diminished by various 
three-dimensional arrangements of wings and bodies. These examples lead 
to a search for some general statements or criteria regarding the wave 
drag of such three-dimensional arrangements. 

In order to put the question in a definite form it will be assumed 
that the airfoils and bodies are contained in the interior of a definite 
three-dimensional region R. The total lift L and the volume v are 
assumed to be given. It is supposed that the wave drag D depends some- 
how on the distribution of the. lift and the volume throughout R and that 
with distributions of a certain family (called "optimum” ones) the drag 
will have a minimum value. It is desired to find the optimum distribu- 
tions, or the conditions determining them, and the value of the minimum 
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drag. Problems of this type have been considered by E, W. Graham and 
his colleagues who give, for example, the optimum distribution of lift 
within a spherical region. 

If the region R is restricted to the plan form S of a planar wing 
then problems of a type previously discussed by the present writer are 
obtained (refs. 5 6 ). In connection with the latter problems it was 

found that all distributions of lift or volume satisfying the given 
requirements could be characterized by relatively simple conditions. The 
present paper describes briefly the extension of these conditions to three 
dimensional regions and the additional conditions required. 

As is usual in linearized -flow problems it will be assumed that the 
disturbance field of the airfoils and bodies can be produced by the action 
of a distribution of sources and ^’lifting elements" or horseshoe vortices- 
One of the difficulties associated with these problems is the determina- 
tion of the actual geometric shapes produced by the distribution of singu- 
larities. In the present analysis the relation between the body shapes 
and the singularities is not known nor determined in detail. For slender 
bodies or thin airfoils closed within the region R it can be assumed 
that the total volume is proportional to the first moment of the source 
distribution with respect to a plane perpendicular to the flight direc- 
tion, whereas the total lift is proportional to the total strength of 
the lifting elements . 

Suppose a region R together with a distribution of singularities 
such as sources or lifting vortices is given. (See fig* !•) Then by 
Hayes ^ theorem (ref. 7)^ the drag will be unchanged by a reversal of the 
whole system. (The geometry of the flow, Including that of the airfoils 
and bodies, will be changed by the reversal but the total lift and the 
total volume will not.) Then the drag may be computed by means of a fic- 
titious "combined disturbance field" obtained by superimposing the dis- 
turbances in the forward and the reversed motion. The perturbation veloc- 
ities in this combined field may be denoted by 

2 u = U|. + Up 
2 v = v^ + Vp 
2 w = W|. + Wp 


An arrangement of sources or lifting elements or their combination which 
yields the minimum drag is then characterized by the conditions 
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w = Constant 


= Constant 
ox 


( 1 ) 


throughout R. 

If conditions (l) are satisfied, then the integrated drag of the 
whole system will he given simply hy 




ax 


The first term on the right-hand side of this expression will he recog- 
nized as the drag arising from a rearward inclination of the lift vector, 
whereas the second term is simply the product of the volume and the con- 
stant gradient of pressure in the combined flow field. 

These conditions may he verified hy making use of a "mutual drag 
relation," essentially similar to the well-known Ursell-Ward reciprocal 
relation, which connects the drag of ary two Interfering distributions 
of singularities in the combined flow field. According to this relation 
the drag of distribution A caused hy the interference of a second dis- 
tribution B is equal to the drag added to B hy the interference of A. 
Now let A he a distribution within R^ satisfying conditions (l). For 

B select a distribution having zero total lift and zero total volume. If 
Rp is contained within R^, then the addition of B will amoruat simply 

to a redistribution, without changing the total lift L or the volume v, 
of A. The drag of A + B may then be written in shorthand notation 

D(A + b) = + Dgg 

Then, since by the mutual drag relation is equal to Dg^, this 

equation may be written as 

D(A + B) = + 2Dg^ + Dgg 

Here Dg^ is the drag of B in the combined disturbance field of A. 

Since w = Constant, v = 0, and = Constant in Ra, this interfer- 

ox 


ence drag may be written as 
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However, since and are zero, vanishes and the added dra^ 

is that of distribution B alone or Dgg. Now the drag of a system alone, 

that is, without interference, cannot be negative; hence, D(a + B) can- 
not be less than D(a) under conditions (l). 


On the other hand, suppose, for example, that the sidewash v^ were 

not zero. A distribution of lateral forces could then be found which 
would result in a negative interference drag, dominating the quadratic 
term D-gg, so that the total drag could be reduced. Hence, if the drag 

of distribution A actually is a minimum value, conditions (l) must be 
complied with. 


Since w 




M ba. 

— and — = — it can be seen that condi- 






tions (1) do not agree with the linearized flow equation 


fl - ^ = 0 

V /^xx ^yy ^zz 


in general, but only if — = 0 . Since — is proportional to the drag 

bx ox 

per unit volume, one concludes that the drag cannot be minimized in an 
absolute sense unless the drag associated with the volume of the system 
is zero (or unless the distribution of singularities is continuous 
throughout R) . Examples such as the Busemann biplane satisfy the former 

condition, namely, — = 0. 

bx 

As Graham et al. have pointed out, distributions of the sort being 
considered here are not unique, since other solutions such" as those shown 
in figure 2 may be added to them without changing the lift or the drag. 

It is interesting to note that conditions analogous to the condi- 
tions w = Constant and v = 0 were found by Monk in connection with 
the vortex drag of lifting systems at subsonic speeds. In that problem 
the conditions apply to the two-dimensional motion associated with the 
trace of the wing system in the Trefftz plane. If the idea of super- 
imposed flow fields is utilized in the subsonic problem, one finds that 
the cylindrical flow associated with the Trefftz plane extends along the 
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whole flight path and hence includes the region Conditions (l) thus 

apply at both subsonic and supersonic speeds^ but are unnecessarily 
restrictive at subsonic speeds * 

Munk’s conditions of constant downwash and zero sidewash were used 
by Hemke (ref, 8) to determine the effectiveness of end plates in reducing 
the vortex drag of a wing at low speeds . It will be interesting to see 
how the condition v = 0 might be used to determine an optimum setting 
and camber for such a surface under more general conditions . This appli- 
cation is illustrated in figure 3 for an end plate on the tip of a wing* 

With the wing in forward motion^ the lateral velocity v^ at the 
surface of the end plate is simply the lateral slope of the fin surface 
times the stream velocity. The condition v = 0 implies that v^ - 
and this condition is obviously satisfied by keeping the geometry of the 
fin fixed when the flow is reversed. At the same time^ however^ recall 
that the distribution of lift and lateral force must be kept the same in 
forward and reversed flow. Hence ^ in order to achieve the minimum drag 
one must find the particular camber and setting of the fin that will yield 
the same distribution of lateral force for either direction of motion. 

At first it seems impossible to satisfy such a requirement since, for 
example, the direction of lift of an inclined surface is usually reversed 
by reversing the direction of flow. However, the form of the adjacent 
wing surface must, in general, change with reversal, since w ^ 0 and 
since the lift distribution on the wing must remain unchanged. Then it 
is evident that the conditions might be satisfied if the pressures on the 
fin surface were dominated by the wing pressures through interference. 

It must be admitted that the considerations have thus far been rather 
abstract, A more concrete result would yield the actual magnitudes of the 
minimum drag associated with various regions. Such results for distribu- 
tions of lift in spherical and ellipsoidal regions have been given in 
reference 4, A somewhat more general result, applicable to arbitrary 
regions R, can be obtained if merely a lower bound for the wave drag is 
sought rather than the actual minimum value . Since this lower bound coin- 
cides with the minimum value in the examples found thus far, it may be 
taken as an approximation to the actual drag in many cases . 

To obtain such a lower bound, we may use Hayes ^ formula (ref, 7) 
or the formula of Lomax (ref, 9)> which expresses the drag more directly 
in terms of areas and pressures intercepted by characteristic planes. 

By utilizing Hayes ^ method of equivalent positions (ref, 7) or the present 
writer *s method of superimposing plane waves (ref. 6), one can construct, 
at each angle 0, three equivalent linear distributions, namely, a volume 
distribution, a lift distribution, and a' side-force distribution. By a 
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harmonic analysis (ref. 10) it is possible to show that the drag associated 
with the leading term in the expansion of the lift distribution i(x), pro- 
port ioiaal to the total lift L, cannot be diminished by interference. The 
possibility^ already known, that the drag associated with the volume can 
be eliminated appears in this analysis. Hence, for the lower boimd the 
value given by the first term in the expansion of the lift distribution is 
used. This step amounts to the assumption that each "lifting line" obtained 
by integrating the spatial lift distribution over the intersecting Mach 
planes is elliptically loaded. For a single elliptically loaded lifting 
line parallel to the flight direction, the wave drag is 


%AVE “ 




( 2 ) 


where I is the length of the line. For the whole region R the fol- 
lowing is obtained 


%AVE = 




where 


(3) 


pi 2 jc o 

_ 1 / sin ‘^9 d 9 /l 

and Z( 0 ) is the projected length of R as defined in figure with 

p = 7 m2 - 1. 

It will be evident from equation (5) that the wave drag depends 
Inversely on the sq-uare of an average projected length of the airfoil 
system - just as the vortex drag depends inversely on the square of the 
span. However, because of the weighting factor sin ^9 the lateral 
dimensions of R are relatively imimportant compared to the dimension, 
or length, along the flight direction. Figure 5 shows the magnitude 
of the error made by using the actual length I and equation (2) for 
the wave drag of several lifting svirfaces. 


Generally speaking, the losses associated with the production of a 
given force in a frictionless fluid are diminished by Increasing the area 
involved in the production of the force and diminishing the pressure. 

Thus the wave drag is diminished by making the "area" as large as 
possible. The vortex drag is diminished by making the square of the span 
as large as possible. On the other hand, to diminish the friction drag 
the actual area S of the wing system must be made as small as possible. 
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At subsonic speeds the conditions are satisfied by making b^ large com- 
pared with S or using a wing of high aspect ratio. It is a matter of 
ordinary observation to see that this condition determines the rather 
special form of subsonic aircraft. At supersonic speeds^ a large value 


of the "longitudinal aspect ratio" 



is needed in addition. 


At subsonic speeds, the elliptically loaded lifting line achieves 
the minimum value of the pressure drag for the whole area covered by the 
wake of the lifting line. At supersonic speeds such a lifting line 
develops, according to linear theory, an infinite drag. However, if the 
line is yawed behind the Mach angle the drag is finite and is actually 
the smallest value obtainable by any distribution within the region of 
the parallelogram ABCD shown in figure 6. Such an oblique lifting line 


maximizes both — and — 


simultaneously . 


At moderate supersonic Mach 


numbers, the results obtained with a V-shaped lifting line - approximating 
a swept wing - are nearly as good. 


When a wing is made narrower so as to approach a "lifting line" while 
maintaining a fixed total lift, the lifting pressure must increase. Even- 
tually the* pressure, or the lift coefficient, will exceed the limitation 
imposed by the small-disturbance theory, or flow separation will occur. 
Beyond this point increases of aspect ratio either laterally or longi- 
tudinally will not necessarily increase the lift-drag ratio. 


Ames. Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Nov. 1, 1955 
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CONDITIONS FOR MINIMUM DRAG 
DISTRIBUTIONS OF LIFT AND VOLUME IN REGION R 



Figure 1 


DISTRIBUTIONS OF LIFT AND VOLUME WITH 
SELF-CONTAINED WAVE SYSTEMS 



UNIFORMLY LOADED 
DISK +L 





ELLIPTICALLY LOADED 
LIFTING LINE -L 


Figure 2 
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USE OF CONDITION v =0 TO DETERMINE OPTIMUM SETTING 
OF VERTICAL FIN ON WING TIP 



APf = Ap^ ; LATERAL FORCE DISTRIBUTION UNCHANGED 
Vf = -V^ ; FIN GEOMETRY UNCHANGED 


Figure 5 


LOWER BOUND FOR WAVE DRAG ASSOCIATED 
WITH THE REGION R AND THE LIFT L 



Figure L 
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APPROXIMATE EXPRESSION FOR WAVE DRAG OF 
LIFTING SURFACE 



m2-i 

Dwave ■ K - 


_Li_ 

TTql^ 



Figure 5 


IDEAL DISTRIBUTION OF LIFT 
FOR PARALLELOGRAM ABCD 
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